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Abstract 

The paper investigates the sensitivity of the inverse problem of recovering the velocity 
field in a bounded domain from the boundary dynamic Dirichlet-to-Neumann map (DDtN) 
for the wave equation. Three main results are obtained: (1) assuming that two velocity 
fields are non-trapping and are equal to a constant near the boundary, it is shown that 
the two induced scattering relations must be identical if their corresponding DDtN maps 
are sufficiently close; (2) a geodesic X-ray transform operator with matrix-valued weight is 
introduced by linearizing the operator which associates each velocity field with its induced 
Hamiltonian fiow. A selected set of geodesies whose conormal bundle can cover the cotangent 
space at an interior point is used to recover the singularity of the X-ray transformed function 
at the point; a local stability estimate is established for this case. Although fold caustics 
are allowed along these geodesies, it is required that these caustics contribute to a smoother 
term in the transform than the point itself. The existence of such a set of geodesies is 
guaranteed under some natural assumptions in dimension greater than or equal to three by 
the classification result on caustics and regularity theory of Fourier Integral Operators. The 
interior point with the above required set of geodesies is called "fold-regular" ; (3) assuming 
that a background velocity field with every interior point fold-regular is fixed and another 
velocity field is sufficiently close to it and satisfies a certain orthogonality condition, it is 
shown that if the two corresponding DDtN maps are sufficiently close then they must be 
equal. 



1 Introduction 



This paper is concerned with the sensitivity (or stability) of the inverse problem of recovering 
the velocity field in a domain from the boundary dynamic Dirichlet-to-Neumann map (DDtN) 
in the wave equation. Let Q he a bounded strictly convex smooth domain in M"^, d > 2, with 
boundary T. Let c(x) be a velocity field in which characterizes the wave speed in the medium 
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and let T be a sufficiently large positive number. We consider the following wave equation 
system: 

-^utt-/^u = 0, gM-^x (0,T) (1) 

u{0,x) =ut{0,x) = 0, xen, (2) 
u{x,t) = f{x,t), gTx (0,T). (3) 

For each / G Hq{[0,T] x F), it is known that (see for instance |16j ) there exists an unique 
solution u E C^{0,T;L'^{n))f]C{0,T;H^{n)), and furthermore |^ S L'^{[0,T] x T), where u is 
the unit outward normal to the boundary. The DDtN map Ac is defined by 

^cUJ := ^l[o,T]xr- 

The inverse problem is to recover the velocity function c from the DDtN map A^. The uniqueness 
of the inverse problem is solved by the boundary control method ffi'st introduced by Belishev 
in [7j. The method can also be used to solve the uniqueness for more general problems, for 
instance, the anisotropic medium case. We refer to [8], [10], [9], [15] and the references therein 
for more discussions. 

We are interested in the sensitivity question for the above inverse problem. Namely, we want 
to investigate how sensitive or stable is it to recover the velocity field from the DDtN map and 
characterize how a small change in the DDtN map affects the recovered velocity field. 

The inverse problem of recovering velocity field is closely related to the inverse kinematic 
problem in geophysics, we refer to [21] for more discussions on the topic. It also can be viewed 
as a special case of the inverse problem of recovering a Riemannian metric on a Riemannian 
manifold. Indeed, it corresponds to the case when the metrics are restricted to the class of those 
which are conformal to the Euclidean one. The inverse problem of recovering a Riemannian 
metric has been extensively studied in the literature. The uniqueness is proved by Belishev 
and Kurylev in [lOj by using the boundary control method. However, as pointed out in [26], 
their approach is unlikely to give a stability estimate since it uses in an essential way an unique 
continuation property of the wave equation. 

The first stability result on the determination of the metric from the DDtN map was given by 
Stefanov and Uhlmann in [23], where they proved conditional stability of Holder type for metrics 
close enough to the Euclidean one in for ^ 1 in three dimensions. Later, they extended the 
stability result to generic simple metrics, |26j . An important feature of their approach is to first 
derive a stability estimate of recovering the boundary distance function from the DDtN map and 
then apply existing results from the boundary rigidity problem in geometry. Their approach 
was extended by Montalto in [M] to study the more general problem of determine a metric, 
a co-vector and a potential simultaneously from the DDtN map, and a similar Holder type 
conditional stability result was obtained. The stability of the inverse problem of determining 
the conformal factor to a fixed simple metric was studied by Bellassoued and Ferreira in |llj . 
They proved the Holder type conditional stability result for the case when the conformal factors 
are close to one. We comment that the result in [11] holds for all simple metrics. For other 
stability results on the related problems, we refer to the references in [34]. 

We emphasize that all of the above stability results deal with the case when the metrics 
are simple. To our best knowledge, no stability result is available in the general case when 



2 



the metrics are not simple. This paper is devoted to the study of the general case when the 
metric induced by the velocity field is not simple. To avoid technical complications due to the 
boundary, we restrict our study to situation when the velocity fields are equal to one near the 
boundary. From this point of view, our results can be regarded as interior estimates. We refer 
to |26j . |28j and the references therein for useful boundary estimates. 

We now give a brief account of the approach and results in the paper. We first derive a 
sensitivity result of recovering the scattering relation from the DDtN map. Our result shows that 
two scattering relations must be identical if the two corresponding DDtN maps are sufficiently 
close in some suitable norm. Equivalently, any arbitrarily small change in the scattering relation 
can imply a certain change in the DDtN map. To our best knowledge, this seems to be the 
first sensitivity result for the problem in the non-simple metric case. Moreover, our result 
is fundamentally different from those in the literature where Lipschitz, Holder or logarithmic 
estimates are derived , see for examples [6], [I], |3l], [23], |2l], [25], [26], [M] and [H]. This 
is the reason the term "sensitivity analysis" is used instead of "stability estimates" throughout 
the paper. We remark that when the geometry induced by the velocity field is simple, the 
scattering relation is equivalent to the boundary distance function. In that case, a Holder type 
interior stability estimate for recovering the boundary distance function from DDtN map has 
been established in [26]. Compared to the Holder type result, our result is much stronger. Our 
approach is based on Gaussian beam solutions to the wave equation, which are capable of dealing 
with caustics, major obstacles to the construction of classic geometric-optics solutions. We refer 
to [20] and [15] for more discussions on Gaussian beams and its applications. 

We observe that for any velocity filed c, the induced Hamiltonian flow T-L^. when restricted to 
the unit cosphere bundle S*M.'^ determines the scattering relation ©c- We linearize the operator 
which maps c to T-L*^\gtj^d and obtain a geodesic X-ray transform operator 3c with matrix- valued 
weight. Note that the scattering relation (or Hamiltonian flow) is the natural object to study 
when the metric is not simple. It is related to the lens rigidity problem in geometry. We refer 
to [30] for more discussions on the topic. The boundary distance function (global or local) has 
received extensive attention in the literature, whose has revealed many important results in the 
case of simple metrics and regular metrics, see for instance [25], [30]. However, this approach 
is unlikely to work in the case of general non-simple metrics. In this paper, we attempt to 
overcome the difficulty by analyzing the scattering relation (or Hamiltonian flow). 

We study the inverse problem of recovering a vector-valued function / from its weighted 
geodesic transform 3cf- For a fixed interior point x, we use a carefully selected set of geodesies 
whose conormal bundle can cover the cotangent space T*W^ to recover the singularity of / at 
X. We allow fold caustics along these geodesies, but require that these caustics contribute to a 
smoother term in the transform than x itself. It is still an open problem to show that such a set 
of geodesies exists generically for a general velocity field with caustics. But we draw evidence 
from the classification result on caustics and regularity theory of Fourier Integral Operators 
(FIOs) to show that it is the case under some natural assumptions in the dimensions equal or 
greater than three. We call the interior point with the above set of geodesies "fold-regular" . A 
local stability estimate is derived near a fold-regular point. 

Finally, we combine the stability result on the X-ray transform and the sensitivity result 
on recovering the scattering relation from the DDtN map to obtain a sensitivity result for the 
inverse problem. 

The paper is organized as follows. In Section 2, we present some preliminaries. In Section 3, 
we introduce the main results in the paper. Section 4 is devoted to the construction of Gaussian 
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beam solutions to the wave equation. The Gaussian beam solutions are used in Section 5 to 
prove the sensitivity result of determining scattering relation from DDtN map. In Section 6, 
we discuss the concept "fold-regular" and prove a local stability estimate for the geodesic X-ray 
transform 3c- In Section 7, we prove the sensitivity result of recovering the velocity field from 
the DDtN map. 

Throughout, we use the following conventions: 

1. Let / and g be two elements in a Hilbert space, then {f,g) stands for their inner product; 

2. Let Ml and M2 be two matrices (including vectors which can be regarded as single column 
or single row matrices), then the product of Mi and M2 is denoted by Mi ■M2. Sometimes, 
the dot is omitted for simplicity; 

3. Let M be a matrix, then M^ stands for its transpose. The same applies when M is a linear 
operator. If M is real and symmetric and C a real number, then M > C means that the 
matrix M — C ■ Id is symmetric and positive definite. If M is a complex matrix, then we 
use for its real part and for its imaginary part; 

4. Let U and V be two open set in a metric space, then U <^V means that the closure of U, 
denoted by U is compact and is a subset of V; 

5. Let Ci and C2 be two positive numbers, then Ci < C2 means that Ci < C ■ C2 for some 
constant C > independent of Ci and C2. 

2 Preliminaries 

In this section, we introduce some notations and definitions. Let J7 be a strictly convex smooth 
domain in with boundary T. Let c be a smooth velocity field defined in Q which is equal to 
one near the boundary. Then c has natural extension to M'^. Throughout the paper, we always 
use the natural coordinate system of the cotangent bundle T*M"' in which we write (x,^) for 
the CO- vector ^jdx^ in T*M.'^. For ease of notation, we also use ^ for the co- vector ^jdx^ . The 
meaning of ^ should be clear from the context. The velocity field c introduces a Hamiltonian 
function Hc{x,^) = ^c^(a;)|^p to T*W^. It also defines a norm to each cotangent space T*W^ by 



Throughout the paper, | • | stands for the usual Euclidean norm in M*^, while | • |c stands 
for the norm in T*W^ induced by the function c. When there is no other velocity field in the 
context, we drop the subscript c and write || • || instead. 

Denote the corresponding Hamiltonian flow by "H*, i.e. for each (xq, Co) £ T*W^, T-L^xq, ^q) = 
{x{t,xo,^o),^{t,xo,(,o)) solves the following equations: 



|e|c = c(x)|e|, for CET^M'^. 



dHc 



(? ■ x{0) = Xq, 



(4) 



X 



dHc 



dx 



(5) 
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We call {x{-,xo,Co),(,{-,xo,^o)) the bicharateristic curve emanating from (xo,Co) and x{-,xo,S,o) 
the geodesic. By the assumptions on c, the flow 7^* is defined for all t G M. Note that the flow 

is also well-defined on the cosphere bundle S*W^ = {(x,^) : x G M, |^|c = 1}. 

We say that a velocity field c is non-trapping in ft for time T > if the following condition 
is satisfied: 

n'^{s*n)f]s*n = (i>. (6) 

Denote 

S;r = {{x,O:xeT,\^\, = l,{^,u{x))>0}; 
SIT = {{x,0:xeT,\^\, = l,{^,u{x)) <0}. 

Assume that the velocity field c is non-trapping in J7 for time T; we now define the scattering 
relation &c : S'l^T — )• 5!^r. For each (xo,^o) £ let /(xo,Co) be the first moment that the 

geodesic x(-,xo,^o) hits the boundary T. Define 

For future reference, we define /_ : S*Q — )• (— oo,0] by letting /_(x,^) be the first negative 
moment that the bicharacteristic curve ?^*(rc,^) hits the boundary S^F and r : S*Q — )■ SZ^T by 

We remark that /-(•) and r(-) are well-defined by the assumption ([6]). 

We now introduce the class of admissible velocity fields that are considered in the paper. 

Definition 2.1. Let Mq, and T be positive numbers. A velocity field c is said to belong to 
the admissible class 2t(Mo, eo, f^, T) if and only if the following three conditions are satisfied: 

1. ce C^{R'^), < ^ < c < Mo, and ||c||c3(iRd) < Mq; 

2. the support o/c — 1 is contained in the set $7^,, =: {x G : dist{x,T) < cq}; 

3. the Hamiltonian He is non-trapping in for time T . 

By Condition 2 above, we can find two small positive constants e* and ei, both depending 
on eo, such that for any (xo,^o) e 'S'*r, if {'H*(xo,^o) : t G (0, /(xq, ^o))} fl 5'*57,o / then 

(Co,i^(xo)) < -e*, (7) 
{ii,y{xi)) > e*, (8) 

/(xo,eo) > ei, (9) 

where (xi,^i) = &c{xo,^o)- 

Finally, we remark that we set up the discussion in the paper in the cotangent space T*W^. 
But one can also set up the discussion in the tangent space TM'^, see for instance [22], [30]. The 
equivalence of the two setups can be seen from the procedure of "raising and lowing indices" in 
Riemannian geometry. We choose the cotangent setup mainly because the following three rea- 
sons. First, it is more natural to the construction of Gaussian beams. Second, the classification 
result of singular Lagrangian maps is more complete than that of singular exponential maps in 
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the literature, though these two problems are equivalent in Riemannian manifold. Finally, it is 
more natural to study caustics in the cotangent space. 



3 Statement of the main results 

3.1 Sensitivity of recovering the scattering relation from the DDtN map 

It is known that the DDtN map Ac determines the scattering relation &c uniquely We 
show that the following sensitivity result of recovering the scattering relation from the DDtN 
map holds. The proof is given in Section [5l 

Theorem 3.1. Let c and c be two velocity fields in the class 2l(eo, ^2, Mq, T). Then there exists 
a constant 5 > such that 

if l|Ac - Ac||/^l[o,3ei/4]xr^L2([0,T+ei]xr) < ^■ 

Remark 3.1. The same result holds when the velocity fields are replaced by symmetric positive 
definite matrices. 



3.2 Linearization of the operator which maps velocity fields to Hamiltonian 
flows 

We begin with the following observation. 

Lemma 3.1. Let c and c be two velocity fields in the class 2t(eo, ri, Mq, T), then ©c = ©c if cLnd 
only ifnl\s''_r = ^I'lsir- 

The above lemma shows the equivalence of the Hamiltonian flow and the scattering relation. 
The next lemma shows that 'H* satisfies an equivalent ordinary differential equation (ODE) 
system in S*W^. 

Lemma 3.2. Let (xcCo) e S*W^ = {{x,i) € M.^'^ : c{x)\^\ = 1}, and let = ^^(xc^o), 

then {x{t),S,(t)) satisfies the following ODE system 

^ = (10) 
i = b{x). (11) 

where b{x) = -iVlnc^. Conversely, if {x{t),^{t)) G S*R'^ satisfies the ODE system / f7^) -/f77]]. 
then {x{t),m) =K{^o,Co)- 

We next linearize the operator which maps each velocity field to its induced Hamiltonian flow 
restricted to the cosphere bundle. Let c be a fixed smooth background velocity field. Denote the 
perturbed velocity field and Hamiltonian flow at time T as = + and = TiJ + 
respectively. Denote also that 56 = — ^V(lnc^ — Inc^) and 



For each {xo,^o) E 5ir, let <I>(t, xo,^o) be the solution of the followmg ODE system 

^{t) = -^{t)A{'Hi), $(0) = Id. 

By the results in Appendix 18. H we have 

x-ijT 

5n^='-^{6b)+r{5b), 

where 



xn/T f-T / n \ 

HSb)M) = l ^"HT,xo,eo)-<^(.,xo,eo)( j ds (12) 



6b 



and < C||(56||^i for some constant C > depending only on ||c||c3( 

Formula (112^ motivates us to define the following geodesic X-ray transform operator 

Mf)i^o,Co)= r Hs,xo,Co)fixis,xo,^o))ds, f G£'in,R^''). (13) 
Jo 

Then ^(56)(rEo,^o) = ^-HT,xo,^o) ■ Uf){xo,^o) with 

^ = ( iV(lnc2°-lnc2) )• ^^^^ 

We associate each (x,^) G a matrix <I>(x,^). Let (^0,^0) = = '^^"''^'^^(x, ^). We 

then define 

<^{x,0 = H-l-{x,0,r{x,0)- 
It is clear that the following identity holds 

^K{xo,^o)) = Hs,xo,^o) 

for all s G M+ such that ^c(^O)Co) £ 5'*r2. We can rewrite the X-ray transform operator 3c in 
the following standard form 

Jo 

HK{xo,Co))f{TT{Kixo, ^0))) ds. (15) 

Remark 3.2. Formula /il5]) is derived in the coordinate ofT*M.'^. Hence it may not be geomet- 
rically invariant. 

Lemma 3.3. Assume that &c = &c, let f be defined as in then 

Pcfh^ < Wffc^ny 
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3.3 Fold-regular points and local stability for geodesic X-ray transform 

We consider the stability estimate of the operator 3c- For simphcity, we drop the subscript c. 
Define /3 : T*R'^\{{x,0) : x e R'^} ^ S*R'^ by 

/3(-,0=(x.||j). 

Let vr : T*M'^ — ^ be the natural projection onto the base space. We define (/> : T*]R'^ — ^ R*^ 

by 

ct){x,0=Tron'=\x,0, {x,0(^T*R''. 

We remark that <j) defined above is equivalent to the exponential map in Riemannian manifold. 
The following result about the normal operator 01 = 3^3 is well-known. 

Lemma 3.4. The normal operator : M^*^) — )• L^(il, M^'^) is bounded and has the following 

representation 

mf{x)= [ W{x,Of{(l){x,0)d(T,{0, /GL2(J7,M2rf) (16) 
JT^n 

where dax denotes the measure in the space T*M'^ induced by the velocity field c, i.e. dcjx{^) = 
c{xYd^, and W is defined as 

^(^'0 = ^pr^*^^ ° • ^ ° ^ ° + ° -0 • ^ ° /3 ° -0}- (17) 

Proof. See [33] or [2i]. 

We see from ()16p that the local property of the normal operator 91 restricted to a small 
neighborhood of x G i7 is determined by the lagrangian map (t){x, •) : T*W^ — )• W^. When the 
map is a diffeomorphism, it is known that the operator ^ near x is a pseudo-differential operator 
(^'DO). However, in general case, the map may not be a diffeomorphism and may have singular 
points which are called caustic vectors. The value of the map at caustic vectors are called 
caustics. When caustics occur, the Schwartz kernel of the operator 91 has two singularities, 
one is from the diagonal which contributes to a ^'DO and the other is from the caustics 
which contributes to a singular integral operator '^2- The property of depends on the type 
of caustics. The case for fold caustics is investigated in [33], where it is shown that fold caustics 
contribute a Fourier Integral Operator (FIO) to 0l2- Little is known for caustics of other type. 
Here we recall the following definition of fold caustics. 

Definition 3.1. Let f : — >■ be a germ of C°° map at xq, then xq is said to be a fold 
vector and /(xq) a fold caustic if the following two conditions are satisfied: 

1. the rank of df at xq equals to n — 1 and det df vanishes of order 1 at xq; 

2. the kernel of the matrix df{xo) is transversal to the manifold {x : det df{x) = 0} at xq. 

We now introduce the following concept of "operator germ" to characterize the contribution 
of an infinitesimal neighborhood of a caustic or a regular point to the normal operator 9T. 
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Definition 3.2. For each ^ G r*M \0, the operator germ 91^ is defined to be the equivalent class 
of operators in the following form 

'yidiy)= [ W{y,v)fWy,v))x{y,v)day{v). (18) 

where x is a smooth function supported in a small neighborhood of (x,^) in M^"^. Two operators 
with xi dnd X2 said to be equivalent if there exists a neighborhood i?(x, ^) of {x, ^) such that 
Xi=X2- X3 for some xs G C^{B{x,^)) with X3{x,0 / 0. 

The operator germ is said to has certain property if there exists a neighborhood B{x,^) 
of {x, ^) in T*W^ such that the property holds for all operators of the form il8\) with x ^ 
C^{B{x,0). 

Properties of the above defined operator germ will be given in Section 16.11 
We note from the preceding discussion that it is complicated to analyze the full operator ^ 
which contains information from all geodesies. However, for a given interior point x, to recover / 
or the singularity of / at x from its geodesic transform, we need only to select a set of geodesies 
whose conormal bundle can cover the cotangent space r*R'^. Caustics may be allowed along 
these geodesies as long as they are of the simplest type, i.e fold type so that we can analyze 
their contributions. This idea can be carried out by introducing a cut-off function for the set 
of geodesies as we do now. We remark that this idea is motivated by the work |28j . For any 
a G C^{S*_T), we define 

^af{xo,^o) = a{xo,^o) a>(H^(xo,Co))/(vr(H^(xo,Co)))ds (19) 







where (xo,^o) £ 5"* F. Let a" be the unique lift of a to S*Q which is constant along bicharac- 
teristic curves, i.e. a^(x,.^) = a o t(x,^) for (x,^) G S*^l. Then is smooth in S*^l and we 
have 

^af {X0,C0)= {J ■<^){K{xo,Co))f{7T{n',{xo,Co)))ds (20) 



JO 

With the original weight $ being replaced by the new one a" • we similarly can define 
In fact, it is easy to check that is defined as in (I17p with W being replaced by 

Wa{x,0 = j^\aoTol3{x,0\'<^^ol3{x,0-^o(3on{x,0 

° r o /3(x, -C)P$t o /3(x, -C) • o /3 o n-\x, -0- 

It can be shown that with properly chosen a, the analysis of the operator DTq, becomes 
possible and we can recover the singularity of / from 91^/. 

We now give two definitions whose discussions are postponed to Section 6. 

Definition 3.3. A fold vector ^ G T*W^ is called fold-regular if there exists a neighborhood U{x) 
of X such that the operator germ is compact from L'^{Q,^q,'M?'^) to //^([/(x), M^*^) (or from 
//^(fieclR^'^) to H'+^{U{x),R'^'^) for all s G R). 

Definition 3.4. A point x is called fold-regular if there exists a compact subset Z2{x) C S*W^ 
such that the following two conditions are satisfied: 
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1. For each ^ G Z2{x), there exist only singular vectors of fold-regular type along the ray 
{t^ : t E M} for the map 4'{x, ■); 

2. \/ie S'^M^, 3 6* G Z2{x), such that 6* ± ^. 

We remark that Z2{x) parameterizes a subset of geodesies that pass through x and along 
which there exist only fold-regular caustics. 

We now present the main result on the local stability estimate for the geodesic X-ray trans- 
form operator. The proof is given in Section 6. 

Theorem 3.2. Let x* he a fold-regular point, then there exist a cut-off function a € C^{S*LV), 
a neighborhood U{x^) of x^, a compact operator ^2,q from L'^{^l^^^^'^'^) to H'^ {U {x^) and 
a smoothing operator y{ from £' {^l,M?'^) into C°°{U{x^),M.'^'^), such that for any Uq{x^) d U{x^) 
the following holds 

(21) 

for all f G V'{n,^,M.'^'^) and s G M. 

3.4 Sensitivity of recovering the velocity field from the DDtN map 

Definition 3.5. An admissible velocity field c is called fold-regular if all points in il. are fold- 
regular with respect to the Hamiltonian flow T-L\.. 

We have established the following main result on the sensitivity of recovering velocity field 
from DDtN map. For simplicity we only consider the case d = 3, similar results also hold for 
d > 3. The proof is given in Section 7. 

Theorem 3.3. Let c and c he two velocity fields in the class ^{cq,^, Mq,T). Assume that the 
velocity field c is smooth and is fold-regular. Then there exist a finite dimensional subspace 
£ C L^(Qep,M^), and a constant 5 > such that for all c sufficiently close to c in H~{^1) and 
satisfying V{lnc^ - Inc^) _L ||Ac - ^c\\Hl,lo,3ei/4]xr^Lmo,T+ei]xr) < ^ implies that c = c. 

4 Gaussian beam solutions to the wave equation 

Let c be a velocity field in the class 2t(eo, ^2, Mq, T). We construct Gaussian beam solutions to 
the wave equation system ([l|)-([3]) in this section. 

We first construct a Gaussian beam in W^. Following |18j . we define G{x,^) = c(x)|^|. For 
a given (xo,Co) £ SIT, let {x{t),^{t),M{t),a{t)) be the solution to the following ODE system: 



M 



X 



a 





a{to) = A4. 



(22) 
(23) 



The corresponding Gaussian beam with frequency A (A ^ 1) is given as follows 



g{t,x,X) = a{t)e 



,iXT(t,x) 
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where r(t,x) = ^(t) • (x - x{t)) + l{x - x{t))'< M{t){x - x{t)). 

Now, let the beam g impinge on the surface T, we want to construct the reflected beam g~ . 
Without loss of generality, we may assume that the ray x{t) hits T at the point x{ti) = xi. Write 
= Ci- We parameterize F in a neighborhood of xi, say V{xi), by a smooth diffeomorphism 
F : U{xi) —7- V{xi), where U{xi) is a neighborhood of the origin in M*^"^. We require that 
F(0) = xi. With the coordinate x = F{y), we can rewrite functions restricted to the boundary 
r. For example, we rewrite 



g{t,x)=g{t,F{y))=g{t,y), T{t,x) = T{t, F{y)) = f{t,y), for x G 

We next derive formulas for f{t,y) and g{t,y). For this, we need to calculate f(ti,0), 
^(ti,0), ||(ti,0) and A^(ti) =: 4^ 



§7"(*i)0), 1^(^1,0) and M{ti) =: jM-^{ti,0). In fact, by direct calculation, we have 



|(«.0) = -1, |((„0) = (|^(0)).6. 



Moreover, the imaginary part and real part of the matrix ^^^^(^1,0) are given below 
9M(ti) = 



c^{xi)ClQM{h)C, -c^{x,)ClQM{h) ■ f (0) 
-c'(^i)(f (0))t?>M(ti)ei (f (0))t?>M(ti)f (0) 



i?(ti,xi)t9=M(ti)i?(ti,xi] 



iVc2(x)6 -(Vlnc(xi))tf (0) 



mt^) = R{tux,)^^MitMt„x,) + _(a|^(o))tvh;c(x,) ^ ' 0(0)6 J 

where i2(ti,xi) = (c2(xi)6, |f (0)). 
We claim that 

^jM(ti) > 0. 

Indeed, note that the column vectors in the matrix ^(0) are linearly independent and hence 
span the tangent space of the surface F at the point xi. By ([8]), ^ forms a nonzero angle with 
the tangent space and thus is linearly independent with all the column vectors in the matrix 
^(0). Therefore the matrix i?(ti,xi) is invertible, and our claim follows. Using d?]) and ([8]), 
we further have 

9M(ti) > C (24) 

for some C > depending on eo and Mq. 

Now, we have calculated f(ti,0), ^(ti,0), |^(ti,0) and g^^(*i)0)- It follows that 

f{h,y) = f(ti,0) + (^^(ii,0))t(t-ii,y) + (t-ti,y)^^^(tl,0)(t-ti,2/)t 

+0{\{t-t,,y)f) 
a TP 

= ((-1, ^(0)t6), (t - ti,y)) + {t- h, y)M{h){t - h, y)t + 0(|(t - h,y)f). 
dy 
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We proceed to construct the reflected beam g . Write 

5-(t,x,A) = a-(t)e^^^"(*'^) 

with 

r~(t,x) = i-{t) ■ (x - x-{t)) + -(x - x~{t))Ul~{t){x - x~{t)). 

We need to find a~(ti), M~(ti)) such that the g~ + g ^ on the boundary. 

Following [2], we impose the following condition 

dt,yTih,0) = d^^yf'{ti,0), foran|a|<2. (25) 

The above condition with \a\ = gives that f~(ti,0) = f(ti,0) = 0; with \a\ = 1 gives that 

(^(o))tei = (^(o))ter, (26) 

where = ^~(ti). Since the column vectors in (§|-(0))^ spans the tangent space T^-^T, we see 
that the tangential component of and are equal. Besides, note that = = 1. Thus, 

r(ii) = Cr = 6-2(6, i^(a:i))i^(xi). 

Condition ([2^ with |a| = 2 gives that M~{ti) = M{ti). Recall the relation between '^M~{ti) 
and 9M(ti), 5iM~(ti) and KM(ti), we have the following two identities: 

R{ti,xi)'<QM{ti)R{ti,xi) = i?(ti,xi)1"3=M~(ti)i?(ti,xi), 
i?(ti,xi)ts)f?M(ti)i?(ti,xi) = R{ti,xi)^QM~{ti)R{ti,xi) 

/ -lVc2(x)(er-6) \ 

[ 0(o)(6'-6) )■ 

Solving the above equations, we obtain QM^{ti) and 5RM~(ti,) and hence M~(ii). Finally, 
set a~ (ti) = —a{ti). Then all of the four components of (x~ (ti), (ti), (ti), M~ (ti)) are 
constructed. We then solve an ODE system to get (x~(t), ^~ (t), M~(t), (t)) as we did for the 
beam g. This completes the construction for the reflected beam g~ . 

We now present some properties about the constructed beam. The following lemma is crucial 
in the subsequent estimates. We refer to [20] for the proof. 

Lemma 4.1. Both the matrices M{t) and M~{t) are uniformly hounded for t € [0,T + ei]. 
Moreover, there exists C > 0, depending on Mq and eq, such that QM{t) > C and QM~(t) > C 
for all te [0,r + ei]. 

We next introduce two auxiliary beams below 

5*(t,y, A) = a{ti)e'^^*, g:{t,y,\) = a-(ti)e^^^*" , 
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where 

dF 

n = {{-l,^{0)Hi),{t-ti,y)) + {t-ti,y)M{ti,0,0)it-ti,y)\ 
dF 

C = ((-l,^(0)ter),(t-ti,y)> + (t-ti,y)M-(ti,0,0)(i-ti,y)t. 

It is clear that = f ~ and = —g~ . 
Lemma 4.2. 

5(t,y,A) = Ut,yA) + 0{V\) in H\{3ei/4,h + ei/2)xU{xi)), (27) 

g-{t,y,X) = g:{t,y,X) + 0{^) m H\{2,ei/A,h + ei/2) xU{xi)). (28) 

Proof: See Appendix 18.21 

Note that ||5(i, y, A)||/,2((j^_^^/2,ti+€i/2)xC/(a;i)) ~ 1- As a direct consequence of Lemma W?2\ 
we obtain the fohowing norm estimate for the beam g restricted to the boundary T. 

Lemma 4.3. 

\\9{-, ■, A)||L2((ii-ei/2,ti+£i/2)xy{a;i)) ~ 1- (29) 

We now present an H^-norm estimate for g~ + g and an approximation for the Neumann 
data 1^ + 1^ on the boundary. 

Lemma 4.4. 

g-{t,x,\)+9{t,x,\) = 0(Va) «n i7i((3ei/4,ti +ei/2) X y(xi)); (30) 
1^ = 2i\g-{ii,u{xi)) + 0{V\) *n L2((3ei/4,ti + ei/2) xy(xi)](31) 

Proof: See Appendix 18.21 

Now, we are ready to construct Gaussian beam solutions to the initial boundary value 
problem of the wave system ([I])-®. We first choose Xei{t) G Co°(M) such that Xei(i) = 1 
for t G (ei/4,ei/2) and Xe^t) = for t G (-oo, 0) U(3ei/4, oo). Let (x^^o) G SIT and 
(xq,Co) ~ ^ "(^^OiCo) = (a^o — Co)- Let g be the Gaussian beam constructed with the 
initial data 3;(0) = a;Q,^(0) = £,l,M{d) = i ■ Id and a(0) = A3. The beam g is reflected by F at 

= ^c{xo,Co) = T^c'^°'^°\xo,^o) at ti = /(xo,Co) + construct the reflected beam 

g~ by the preceding procedure. Let u be the exact solution to the wave system ([l])-([3]) with 

f{t,x,X) = g{t,x,X) ■ Xeiit). 

Then u = g + g~ + R, where the remaining term R satisfies the following equation system 

VR = -V{g + g'), (t, x) G x (0, ti + ei/2), 

RiO,x,X) = -{g + g^){0,x,X), xen, 

Rt{0,x,X) = -{gt+ gt){0,x,X), x e Q, 

R{t,x,X) = -g{t,x,X){l-Xe^{t))-g~{t,x,X), (t, x) G (0, ti + ei/2) x T. 



13 



Here V stands for the wave operator ^^^-^ du — A. 
Lemma 4.5. 

ll^llL2([o,ti+.i/2]xr) < (-VA 

for some constant C > depending on eo and Mq. 

Proof. We apply Theorem 4.1 in [16] to derive the estimate. Note that the compatibihty 
condition is satisfied on the boundary at time t = 0. It remains to show that the following four 
estimates hold: 

Wna + 9-)\\cm,+e,/2];L^m < A (32) 

\\{g + 9^){0,;X)\\Hi(^n) < A (33) 
\\{gt+gT){0,;X)\\LHn) < A (34) 

\\g{t,X,\){l -Xeiit)) - g~{t,X,\)\\Hi{[0fy+ei/2]xT) ^ (35) 

First, ()32p follows from the standard estimate for Gaussian beams, see for example [5]. We next 
show ()33p . By Lemma |4. 11 there exists a constant C > depending on Mq and eo such that the 
following two inequalities hold 

\g{t,x,X)\ < Al.e-^^-l— 
\g-it,x,X)\ < A3 -e-^^-l^-^'WI'. 

Thus the beam g and g~ are exponentially decaying away from the ray x{t) and x~{t) re- 
spectively. Using this property, it is straightforward to show that 11^(0, •, A)||jji(q) < 1 and 
\\g~{0,-A)\\m(n) ^ 1, whence §S) and §^ follows. 

Now, we show (|35p . We divide the domain (0,ti + ei/2) x P into three parts: 

El = (0,ei/2) X r, S2 = (ei/2,ti-ei/2) xr, S3 = (ti - ei/2, ti + ei/2) x L. 



We show that inequality (j35|) holds on each part. 

For {t,x) G Si, we have 1 — Xtiit) = 0- Consequently, 

g{t,x){l -Xei(t)) -g~{t,x,X) =g~{t,x,X). 
By the exponential decaying property of g~ , we obtain that 

\\g{t,x){l-x,i{t)) - g~{t,x,X)\\Hi(j:,) < Vl- 
For (t, x) € S2, by the exponential decaying property for both g and g~ again, we obtain 

\\g{t,x,X){l - Xeiit)) - g~{t,x,X)\\Hi(T:2) ^ 1- 

Finally, for {t, x) £ S3, note that ti - ^ = 1{xq, (,o) + ^ - ^ > We can apply Lemma [Ol to 
the part x G ^(^i) and the exponential decaying property for both g and g~ to the remaining 
part to conclude that 

\\g{t,x,X){l - Xeiit)) - g~ {t,x, X)\\hi(s,a < 
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This completes the proof of (I35p and hence the lemma. 



5 Proof of Theorem 13.11 

Proof of Theorem D For any (xo,Co) G SIT, let = G^{xo,^o) = 'Hc''"'^°\xo,£,o) and 

(Si, 6) = &c{xo,Co) = 'H't'''^"\^o,^o)- We need to show that (/(xq, Co), a^i, 6) = (K^^o, Co), ^i, li) 
if ||Ac — Acll is sufficiently small. We do this in the following steps. 

Step 1. Let ti = 1{xq,(,q) + ^ and ii = /(xo,Co) + x- Without loss of generality, we may 
assume that ti < ti. Let V{xi) be a neighborhood of xi in F which is parameterized by a smooth 
function F : U{xi) — V{xi) as before. We may assume that xi G V{xi). Let xi = F{5y). We 
construct the initial beam g, the reflected beam g~ , the boundary Dirichlet data /, the solution 
u to the wave equation with velocity field c and remanning term R as in the previous section. 
We similarly construct g, g~, u and R to the system with velocity field c and with boundary 
Dirichlet data f = f- 

Step 2. Denote by I{ti, ei/2) the interval (ti — ei/2, ti+ei/2). Since ti < ii and /(xq, Co) ^ ei, 
we have /(ti, ei/2) C (3ei/4, ti + ei/2) and I{ti, ei/2) C (3ei/4, ti + ei/2). Then we can apply 
(jSip and Lemma 14.51 to obtain 

^ * A N i- du du 
Ag-Ac / = TP 

d{g + g-) d{g + r) , gfi 9^ 
= 2i\ ■ {{^uHxi)) ■ g - (Ci, i^(xi)) • ~g} + 0{^) 

in L\l{h,ei/2)xV{xi)). 
It follows that 

((As - Ac)/, <7>i2(j(t^,,^/2)xy(:.i)) = 2^A • (6, '^(a^i)) • (5, 5)i2(^(j^,,^/2)xi/(xi)) 

- (iw^ixi)) ■ (5,5>i2(,(t,,,,/2)xF(xi)) +0(\^). 

Note that 

lii^c - -^c)/, g)L^(I{ti,ei/2)xV{xi))\ < ll(Ac - Ac)/||2,2(/(j^^ej/2)xy(xi)) " \\9\\L2{I{ti,ei/2)xV(xi)) 

< \\{Ac - Ac)/||L2((o,T+ei)xr) • \\9\\L^{I{tuei/2)xV{xi)) 

< ||Ac - Ac||^l([o,3ei/4]xr)^L2{[0,T+£i]xr) • ll/llHl([0,3ei/4]xr) 
■\\9\\L2{I{ti,ei/2)xV{xi)) 

< A • ||Ac - Ac||^i([o^3£-^/4]xr)^L2([0,T+ei]xr)- 

Thus the following inequality holds 

\{^l,v{xi)) ■ {g,g)LHlih,ej/2)xVixr))\ " ' {9 , Q) L^{I{t^,e^/2)xV{x^))\ 

C 

^ + • ||Ac - Ac||^i([o^3£-^/4]xr)^L2([o,T+ei]xr) (36) 
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for some constant C > 0. 

Step 3. We now estimate the two terms on the left hand side of the inequality (j36p . First, 
by ([!]) and Lemma l43l we have 

• \{9,9)Lmit,,e,/2)^V{x,))\ ~ 1- (37) 

We next estimate {g, g) L'^[i{ti,t-i_/2)xV{x-i))- the coordinate x = F(y), by Lemma we have 

{g^ g) L'^{i{ti,ei/2)xu{xi)) = {g*^g*)L^{i{ti,ei/2)xu{xi)) "^^("y^)' (3^) 

Here we recall that = o(ti)e*'^'^* and g^ = a(ti)e*'*''^* with 
dF 

n = ((-l,^(0)tei),(t-ti,y)) + (t-ti,y)M(ti,0)(t-ii,y)t; 
dF 

h = ((-1, - ^i^y- ^y)) + - ti,y- dy)M{ti,0){t - ii,y- 6y)l 

By Lemma IS.ll in Appendix 18.31 we have 

\{L,g*)LHlit„e,/2)xUix,))\ <e-'^°^l'^l (39) 

where cq is a positive constant depending only on ||c||c3 + \\c\\c3 and \dz\ = \ti — + + 
\%{Sy)^^i - H (0)"^eiP- It follows from §B) and ([MD that 

\{~g,g)LHiit^,e,/2)xVix,))\ < e-^o^l-^^l +0(-^). (40) 
Step 4. Combining ([MI), (I3ZD and (001), we see that 

g coA|5z| > _ CgllAc - Ac||^i([o_3ei/4]xr)^L2([o,r+ei]xr) - ^3-^ 

for some positive constants Ci, C2 and C3 which are independent of (xq, ^o)- By letting A — )■ 00, 
we conclude that 5z = if 

II A A II ^1 
||Ac - Ac||j^i([0^3^j/4]xr)-5.L2([0,T+ei]xr) < 

Set 5 = g^. Prom 6z = 0it follows that ii = ii, 6y = 0, and ||(0)t|i - ||(0)t^i = 0. It remains 

to show that ^1 = ^i- Indeed, §|-(0)^C^i — §|-(0)^^i = implies that the tangential component 

of ^1 and ^1 are equal. Besides, ||^i|| = ||^i||. These together with yield that ^1 = ^i- This 
completes the proof of the theorem. 
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6 Geodesic X-ray transform with caustics 
6.1 Local properties of the normal operator ^Tt 

In this subsection, we present some results about the local properties of the normal operator 
m{see (HID). 

From now on, we fix G $7. We first decompose locally into two parts based on the 
separation of singularities of its Schwartz kernel. Note that the map (l)ix^,-) : M'' M'^ is a 
diffeomorphism in a neighborhood of the origin. In fact, we can check that ^^^(0) = c{x^)-Id. 
Similar to the proof of existence of uniformly normal neighborhood in Riemannian manifold |17j . 
we can find £2 > and a neighborhood of U{x^) CW^, such that 

4>{x, •)|||5|j<2e2 ^ diffeomorphism for any x G [/(x*). (41) 

Let Xe2 G Cg^(M) be such that x{t) = 1 for |t| < €2 and xit) = for \t\ > 2e2. We then 
define 

= / W{x,0f{4>{x,0)Xe,m\\)da,{0, (42) 

^J{x) = I W{x,i)f{<i^{x,m^-XeM\\))dcT.{i). (43) 

Note that for any / supported in 0, f{4>{x,S,)) = for all ||^|| > T. Thus we have 

^z/(^)= / W{x,C)f{ct>ix,0){l-XeM\\))daAC). 

It is clear that 91/ = 91i/ + Olz/- This gives the promised decomposition of 91. We next 
study Oil and 9I2 separately. 

Lemma 6.1. Ol^ is an elliptic ^DO of order -1 from C^(i/(a;*), M^^) to V {U{x^),M?'^) with 
principle symbol 

c7p(oii)(x,e) = 2^- [ 6{{t0mHx,o)-<^ix,e)da,{e). 

Js^n 

Proof. See [28] or [33] • 

We now proceed to study the operator 9I2 whose property is determined by the Lagrangian 
map (/>(x*, •). We shall study the operator germ ^2,^, for each ^ G T*^M.'^. We first consider the 
case when is not a caustic vector, i.e. ^* is a regular vector. 

Lemma 6.2. Let G 5'*._^M'^ be a regular vector, then there exists a neighborhood U{x^,) of 
and a neighborhood B{x^,S^^) of (x*,^*) such that for any x ^ Co°(-B(a;*, ^*)) the following 
operator 

'yi2,iJix)= [ Wix,C)f{<P{x,m^-XeM\\))-Xix,Od(TAO 

J Tin 

is a smoothing operator from £'{Q.,'M?'^) into C°°{U{x^),'M?'^). 
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Proof. Since G 5'*^M'^ is regular, there exist a neighborhood V{x^,) of in M"' and a 
neighborhood of (x,,,^*) in M^*^ of the form = ^(a;^,) x -Bo(^*) for some open 

set -Bo('^*) in I^'^ such that the map 4>{x, •) is a diffeomorphism between -Bo(C*) a-nd its image for 
all X £ V{x^). We denote the inverse of the map 4>{x, •) by (j)~^{x, •). By a change of coordinate 
^ = cl)~^{x,y) and use some cut-off function, we can write ^2,^, in the following form 

'yi2,iJ{x)= [ K{x,y)f{y)dy,f e£'{n,M.^'') 
Jn 

for some smooth function K in Q x Q. The Lemma follows immediately. 

We next consider the case when is a fold vector. We have the following slightly modified 
result from [33j . 

Lemma 6.3. Let be a fold vector of the map (p{x^:,-). Then there exists a small neigh- 
borhood C/(x*) 0/ X* and a small neighborhood B{x^,£^^) of (x*,.^*) in M^*^ such that for any 
X G the operator 0X2,5, '■ £'{^,'^'^'^) ^ V'{U{xS),M.^'^) defined by 

^2,€J(^)= / W{x,OfiH^,0)a-XeM\\))-x{x,Oda.iO, f€£\n,M.'"') (44) 

is an FIO of order — | whose associated canonical relation is compactly supported in the following 
set 

|(x,^,y,r/); x £ U{x^),y = (I){x,uj),{x,uj) £ B{x*,C*), det (i^(/)(x, w) = 0, 

^ = ^ Cofeer(d^<A(x,a;)).} (45) 

Proof. We sketch a proof here and refer to [33] for detail. We first note that by the fold 
condition, there exists a small neighborhood 5i(x^,,^*) of (x*,^*) in M?'^ such that is the only 
singular vector of the map (/>(x^,, •) along the ray {tS^^ : t G R} in Si(x*,^*). Define 

S = {{x,w) : det d^(f>{x, w) = 0} C M^^, 

S = {{x,y):y = cPix,oj),{x,u;)€S}cR^''. 

By shrinking i?i(x*,,^*) if necessary, we can show that Sf]Bi{x^,^^) is a smooth {2d — 1)- 
dimensional manifold in M^'^, and (/> is a diffeomorphism between S f]Bi{x^,,^^,) and its image. 
Denote 

Si = Sf^Bi{x^,C^), T,i = (p{Si). 

Note that Si is a smooth (2d — l)-dimensional manifold in a neighborhood of (x*,?/*) in x J]. 
Let 7r2 be the projection from x to its second component. By the fold-condition for and 
the fact that the matrix da;,5</'(x*, is surjective, we can show that dTr2 ■ T(^x,,s,,)'^i ~^ Ty^M.'^ 
is surjective. Thus there exists a neighborhood Vi(y*) of y* in R'^ such that ^1(1/*) C 7r2(Si). 
We remark that the surjectivity of the linear map dTr2 implies that the conormal bundle of Si 
belongs to {T*R'^ \ 0) x (r*R'^ \ 0), where stands for the zero section of the normal bundle 

r*R'^. 

Let C/i(x*) be a small neighborhood of x* such that C/i(x*) C vr(S'i) and x ^ C^{Bi{x^,S,*))- 
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Consider the Schwartz kernel of operator 9^2,5. : £'{n,M.'^'^) V'{Ui{x^),R'^'^) defined by 

^2,€J(x)= / Wix,C)ficPix,Cm-XeM\\))-xix,C)da,{C), f € £' in,R^''). 

We can show that it has conormal singularity supported in the set Si. Moreover, the conormal 
bundle A/"* Si is given by 

7V*Si = ^{x,i,y,ri); x £ U{x^),y = (/>(x,a;), {x,u)) E Bi{x^,(,^), 

dtp'' (xj uj^ 1 
^ = -rji , r] G Coker (d(^0(a;,a;)), det d^(j){x,uj) = O.j 

By analyzing the singularity of the Jacobian determinant of d^(f){x,uj)), we can show that the 

Schwartz kernel of ^2,^* belongs to the conormal class x 0, Si, 2^x2^)5 where A^2dx2d 

denotes the vector bundle of matrices from R?'^ to R?'^ over Q.. Especially, when the domain of 
'^2,£_t is restricted to distribution sections supported in Vi{y^), the operator 912,5, ^ of 
order -f from S' {Vi{y^),R?'^) to V'{Ui{x^),R?<^). 

The above define 0^2,5, with x ^ C^{Bi{x^,^^)) requires the domain to be £' {Vi{y^),'M?'^), 
we now show that this condition can be relaxed by further decreasing Indeed, let 

U{x*) be a neighborhood of such that U{x*) d Ui{x.t) and -B(x*,^*) be a neighborhood of 
(x*,^*) such that B{x^,,^^) ^ i?i(x*,^*). By choosing C/(x*) and -B(x*,^*) to be sufficiently 
small, we can assume that the set ^(y*) = { <^(x,^) : x E [/(x^,),^ E i3(x,,,^,,)} is compactly 
supported in Vi{y^). We then choose xv £ C'o°(Fi(x*)) such that xv{x) = 1 for x E y(x*). 
One can check that for any x ^ Co°(-S(^*) ^*))) tlie operator ^2,,^, defined by 

9T2,?J(x)= / VF(x,e)/(<^(x,6)(l-Xe,(||^||))-x(x,0d^^x(0, feS^^^R"") 
Jt*q 

satisfies 9l2,5J(a;) = ^2,iSxv f){x) for ah x E U{x^). Thus 9X2,5. well-defined from £'{n,R^'^) 
to P'(?7(x*), M^'^), and is a FIO of order — | with canonical relation compactly supported in the 
set (j45|) . This completes the proof of the lemma. 

6.2 Singularities of the map (p{x, ■) 

In this subsection, we present some properties about the map (j){x, •) which is equivalent to the 
exponential map in Riemannian manifold. 

By the classification result for Lagrangian maps (see [3] and [3] for detail), there are only a 
finite number of stable and simple singular Lagrangian map germs in dimensions between three 
and five and they are generic. In three dimensions, there are four types: fold, cusp, swallow-tail 
and D4. The others are unstable and can be removed by using arbitrarily small perturbations. 
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We define 

/C(x) = E T*M'^ : the map germ (j){x, •) at ^ is singular}; 

ICi{x) = {(, G T*M.'^ : the map germ (j){x, •) at ^ has singularity of fold type}; 

K^2{x) = {£, £ T*M'^ : the map germ (j){x, •) at ^ has singularity of cusp type}; 
IC3{x) = e T*M.'^ : the map germ (j){x, •) at ^ has simple and stable singularities of types 

other than fold and cusp}; 
ICiix) = G T*R'^ : the map germ (j){x, •) at ^ has singularity which are either not simple or stable}. 

It is clear that /C(x) = Uj=i^j(^)- Denote 

S{x) = {^e S^R'^ : e /C for some r G M}, 
Sj{x) = {£,£ S*R'^ : ri £ ICj for some r £ M}, j = 1, 2, 3, 4. 

We say that the map <j){x, •) is in a general position(or generic) if the map germ <j){x, •) is 
simple and stable at all caustic vectors in IC{x), i.e. /C4(x) = 0. It is possible that the map (j){x, •) 
can be brought to a general position by adding an arbitrarily small perturbation to the velocity 
field c. By the classification result of Lagrangian maps, see for instance [3], the following result 
holds for the set IC{x). 

Propsition 1. Assume that the map (pix,-) is in a general position, then the sets /Ci(x) and 
fC2{x) are smooth manifolds of dimensions d—1 andd—2, respectively. The setlC^{x) is a union 
of smooth manifolds of dimensions not greater than d — 3. Especially, for d = 3, the sets ICi{x), 
IC2{x) and }C-i{x) consists of smooth surfaces, smooth curves and isolated points, respectively. 

In the case when the map (/>(x, •) is not in a general position, it is known that /Ci(x) U fC2{x) |J lC^{x) 
is open and dense in K.{x). 

Note that Sj{x) are the images oiK,j{x) under the map /3 which sends ^ G T*W^ t° M ^ '^^^'^ 
for ^ 7^ 0. We conclude that the following result holds. 

Lemma 6.4. Assume that the map (j){x, •) is in general position, then the sets S2{x) and S^{x) 
are of finite d—2 and d — 3 dimensional Hausdorff measures, respectively. Especially, for d = 3, 
the set S2{x) is a curve (not necessarily smooth) of finite length in S*R^ and Ss{x) consists of 
a finite number of points. 

6.3 Discussions on the concept of Fold-regular 

In this subsection, we discuss the concept "fold-regular". We show that for a general velocity 
field in whose induced metric is not simple, a given point in Q is fold-regular under some 
natural assumptions. 

We begin with the concept "fold-regular vector" . It is still an open problem to find a complete 
characterization for it, i.e. what are the necessary and sufficient conditions for the map germ 
<j){x^:,-) at for to be fold-regular. We have the following partial answer in the form of 
remarks. 

Remark 6.1. In dimension d = 2, the set of fold-regular vectors is generally empty. Indeed, 
for a fold vector the operator germ 9l2,^« is a FIO of order —1, and hence the best estimate 
is that it is bounded from L^(r2ep, M^*^) to i/^(C/(x*), M^*^) for some neighborhood U{x^) of x^,. 
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Remark 6.2. In dimension d > 3, a sufficient condition for a fold vector to be fold-regular 
is that the following condition is satisfied 

dl(j){x*,C*){N^AC*) \ X ■)\t^^s{x,) is of full rank. (46) 

where Nx^{£,=t) denotes the kernel of d^(j){x^,^^) and S{x^) the set of all vectors ^ G T*W^ such 
that det (ig(/)(x*, ^) = 0. Indeed, in that case, it is shown in f3^ that the canonical relation 
associated with the operator germ 912,^, is locally a canonical graph and hence 9T2,^« is bounded 
from L2(Q,Q,M2d) ^ H2 {U{x^),R^'^) for some neighborhood U{x^) of x^. Note that for d > 3, 
H2[U{x^:),R'^'^) is compactly embedded in H^{U{x^:),R'^'^), so 9^2,?. is compact from L'^{Vt^^,R?'^) 
to {U {x^) ,R?'^) and we can conclude that is fold-regular. 

The set of fold-regular vectors Zi[x^) contains more elements than those which satisfy the 
graph condition ^B^. In fact, let C C T*Q. x T*^ be the canonical relation associated with 
the operator germ ^2,5, defined in Lemma \6.S\ We have shown that C is homogeneous and 
C C (T*O\0) X {T*n\0). By the main result in fl^, ^2,^, is bounded from L'^{n^g,R'^'^) to 
if2-3([/(x,),M2d) for some neighborhood C/(x*) of x^,, if the only singularity of the projection 
of C to its first or second component at the point associated with (x*,^*) is fold or cusp. Since 
H2~3{lJ[xjf),R^'^) is compactly embedded in H^{U{x^,),R?'^), we see that'^2,i, is compact from 
L2(Q^g,M2rf) ^Q H^{U{x^),R^'^) and hence is fold-regular. 

We now consider the concept "fold-regular point". We denote 

Zi{x^) = {.^ G 5'*^M'' : V r G M,r^ is either regular or fold-regular for the map(/>(x^,, •)}. 

It is clear that 

Zi{x,) C Z{x,) =: SlX\{S2{x.)\JS:i{x,))\JSi{x,). 

We remark that 2^(x*) characterizes the set of geodesies that pass through x-f and along which 
the map </>(x*, •) only has singularities of fold type. By Morse's index theorem (a fold vector for 
the map (/>(x, •) corresponds to a fold conjugate vector for the exponential map exp^(-)), for each 
^ G 5(x*), there are at most finitely many fold vectors along the geodesic vr o 'H*(x*, ^). Using 
the definition of "fold-regular", we can conclude that 2^i(x*) is open in S'*^M'^. 

Recall that x* is fold-regular if there exists a compact subset Z2 (x* ) d Zi [x^, ) such that the 
following completeness condition is satisfied 

G SIR'^, 30 G Z2{x^), such that 6 ^ i. 

Remark 6.3. A sufficient condition for the completeness of Zi{x^:) is that Zi[x^) contains a 
set 

e^ = {i&siR\iLe} 

for 9 G S* R'^. 

Remark 6.4. // the completeness condition fails for 2i(x*), then there exists 6 G 5';^,IR'^ such 

that 

cS2{x,)\JS3{x,)\JSi{x,). (47) 

Assume that the map (j){x^:, •) is in a general position. By Lemma \6.4\ i/ie sei 52(x*) |J 53(x*) |J 54(x*) 
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is of finite d — 2 -dimensional Hausdorjf measure, so is the set 9^ for each 6 G S'^W^. Thus, we 
conclude that there exists at most a finite number of 6 such that (J^ holds. 

6.4 Proof of Theorem lOl 

We prove Theorem 13.21 in this subsection. The proof can be divided into two major stages: in 
the first stage, we present some preUminaries and construct a cut-off function a S C^(S'ir) 
which selects a complete set of geodesies with only fold-regular caustics, see Lemma l6.5j in the 
second stage, we study the normal operator 9Ta = '^a'^a, see Lemma and \677\ Theorem 13.21 
is then a direct consequence of Lemma 16.61 and 16.71 

We now present some preliminaries that are necessary for the construction of a. Let 
be a fold-regular point with the compact subset ^2(2;*) C 5'*^M'^ in Definition 13.41 Denote 
Ce2,TZ2 = {r^;^ G Z2{x^),r G M and €2 < \r\ < T}. For each G Ce2,T2^2, by Lemma [612] and 
Lemma [Oj there exist a neighborhood U{x^,^^) of and a neighborhood B{x^,S,*) of {x^,S,*) 
such that for any x £ C^{B{x^,^^)) the following operator 

%C./(x)= / W{x,Of{^{x,0)i^-XeM\\))-x{x,Oda,{0 
J Tin 

is compact from H''{Q.^^,M?'^) to C*), M^rf). Let ^0(2;*,^*) be another neighborhood 

of (x*, in M^'^ such that Bq{x^,^^) d ^*). Since Ce2,TZ2 is compact, there exists a finite 

number of .^*'s in C^2,T^2-, say ,^1, ^2vi ^M, such that 

M 

C,2,r-2^2 C IJ So(a;*,Cj)- 
i=i 

We can then find smooth functions xi-, X2,---,XM with suppxj C B{x^,S,j) for each j such that 

M M 
^ Xj{x, = 1 for all (x, G U Bo{x,,Cj). 
j=i j=i 

Denote by ^0 be the greatest connected open symmetric subset in Ujli -^0(2^*; Cj) which 
contains Ce2,T2^2- Here and after, we say that a set B in M^'^ is symmetric if (x,.^) G B implies 
that (x, -C) G B. Define 

A = e IR^"' : k-x,| < e, £2 < lieil <T} 

for each e > 0. It is clear that is compact in M^'^, so is the set Ae\Ao- 

Lemma 6.5. There exist 63 > and a G C^(S'* F) such that the following two conditions are 
satisfied: 

aixo,Co) = 1 /or aZZ (xo,Co) G To/3(Ce2,T-H2(x*)), (48) 

a{xo,Co) = for all {xo,^o) ero 13 {Aes\Ao). (49) 

Proof. Note that both /3 and r are continuous. Since Ce2,T22{x*) and .4e\^o are compact. 
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so are the sets r o /3(Ce2,T^2(a^*)) and r o /3(^^\^o)- We claim that there exists €3 > such that 

r o P{C,,,tZ2{x.)) fl T o /3(A\A) = 
for ah e < £3. Indeed, assume the contrary, then 

for all e > 0. Note that the collection of compact subsets r o j3{Ce2,TZ2{x*)){^T o /3(^^\^o) 
is decreasing with respect to e, so it satisfies the finite intersection property and we can thus 
conclude that 

T o /3(C,,,t22(x*)) fl T o / 0. 

But on the other hand, we can check that 

nro/3(A\A) = T((A\A)n^^'*i^') 

Using the fact that r is injective on 5'*^M'^ and Ce2,T^2 C Aq, we obtain 

r((A\-4o) n SIX) n ^(c..,T^2(x*) n six) = 0- 

Thus, 

T o fi{e,,^TZ2{x,)) f T o = 0. 

This contradiction completes the proof of our claim. 
Now, we have 

r o /3(C,2,t22(x.)) f T o /3(A3\A) = 0. 
By decreasing £3 if necessary, we may assume that 

{x ■.\x — x*| < £3} C 7r(^o)- 

Since both the sets r o j3{Ce2,TZ2{x*)) and r o /3(^e3\-4o) are compact in 5* F, we can find 
a S C(f (5ir) as desired. This concludes the proof of the lemma. 

The construction of a above completes the first stage of the proof of Theorem 13. 2| we are 
now at the second stage. We define the truncated geodesic X-ray transform 3af as in ()19p or 
(j20p . By replacing the weight <^ with the new one • <1>, we obtain OTq, ^i^a and ^2,0 from 
the corresponding formulas of CR, and It is clear that Lemma |6.2| 16.31 still hold with the 
new weight. 

Lemma 6.6. There exist a neighborhood U {x^) ofx^ and a smoothing operator from £'[Q, M?'^) 
into C°°{U{x^,),M.^'^), such that for for any s G M and any neighborhood C/o(x*) 0/ with 
Uo{x^:) ^ [/(x*), the following estimate holds 

ll/ll//=((7o(x.),K2'') ^ ll^l,a/||H''+i((7(x*),M2<') + l|5^/||H=(f1,M2d). (50) 
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Proof. We first show that OTi^q, is an elhptic \I'DO. Indeed, as in Lemma [6. 11 ^i,q is a ^'DO 
of order -1 from C^(i7(x*), M^d) to V'{U{x^),R'^'^) with principle symbol 



By the construction of a, for any ^ G S'*^M'^, we have a^{x^,,9) = 1 for some 9 G 5'*^M"' with 
9 Thus 



in the sense of symmetric positive definite matrix. By continuity, we can find a neighborhood 



C of x^ such that o"p(9Xi^a)(x, > for all x £ U{x^,) and ^ S S'*M'^. Thus we can 

conclude that is an elliptic ^-DO of order -1 from C^([/(x=,), M^^) to P'(f/(x*), M^d). 

Now, let *B be the pseudo-inverse of OTi^a restricted to U{x^:). Then !B is ^'DO of order 1 
and there is a smoothing operator «Hi : £' {U{x^),R^'^) C°°(C/(x*), M^*^) such that 



for all g supported in f7(x*). 

Next, let Uo{x*) be any given neighborhood of x* with Uo{x^) d U{x^). For later conve- 
nience, we write Us{x*) for U{x^). Then there exist two neighborhoods of x*, say C/i(x*) and 
U2{x^) such that C/o(xi) d C/i(x*) d [72(2;*) ^ Us{x^:). We choose three smooth cut-off functions 
Xo, Xi arid X2 such that suppxj C Uj+i and Xjll/j = 1 for J = 0, 1,2. Then both Xo^(l — Xi) 
and xi^i,a(l — X2) are smoothing operators. 

Note that X2/ is compact supported in U3{x^), so we have by ([5T]) that 



Xof = Xo • X2/ = Xo • *B o gii,aX2/ + Xo9^i(X2/) 

= X0«Xl^i,aX2 • / + (X0*B(1 - Xl)51i,aX2 + X0^lX2)/ 

= XO«Xl^i,a/ + (xO*BxiOTx,a(l - X2) + X0«(1 - Xl)9Xi,aX2 + Xo9^iX2)/ 
= X05Sxl^i,a/ + ^/ 

where 9^ = Xo^Xi^i,a(l - X2) + Xo5S(l - Xi)^ i,aX2 + Xo9^iX2- We can check that 9^ is a 
smoothing operator from £'{n,R'^'^) to C°°(C/(x*), M^"'). 





(51) 



X2/ = *Bo9I,,„(x2/) + 5Hl(X2/). 



Thus 



Finally, we conclude that 



\\f\\H=(Uo{x,),R'2'i) 



< 



||X05SXl^l,Q/||Hs((70(i:.),M2d) + \\^f\\H={Uoix^),R'2'i) 
||*Bxi^l,a/||H'>((73(x.),R2'*) + II^/IIh''(C/o(x,),R2'*) 
IIXl ■ ^l,a/ll_H'''+i{!73 + \\m 

||^i,Q/||/f'<+i(;73(x*),R2d) + ||9^/||//s((7o(x.),IR2d). 



< 



< 



< 



This completes the proof of the Lemma. 
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We now study the operator ^2,a- 

Lemma 6.7. There exists a small neighborhood of x^, say ?7(a;*), such that the following de- 
composition holds for the operator 9T2,« : S'{Q.,^^'^) P'(C/(a;*), M^d) 

M 

^2,a = J^0T2,i (52) 

i=i 

where each OTaj is compact from H-'{n,g,R'^'^) to H'+^{U{x^),R'^'^). 
Proof. Recall that ^2,a has the following representation 

'yi2,afix) = [ Wa{x,Of{H^,m^-XeM\\)) ' da^iC), 

J Tin 

where 

Wa{x,0 = j^\aoTol3{x,0\'^^ol3{x,0-^of3onix,0 

By (|49l) and the fact that Aq is symmetric, we see that supp Wa C for all x with |x — x*| < £3. 
Now, let Xj's be as in the first stage. Define 0^2^ : S'{n,R'^'^) V'{U{x^,^j),R'^'^) by 

^2,jf{x) = [ WUx,Of{Hx,Cm-XeM\\))-Xj{x,Oda,{0. 
J Tin 

Let U{x^) = f]'jL^{U{x^,^j))f]{x : \x - x*| < £3}. Then U{x^) is a neighborhood of x* and 
each ^2,j is compact from M^rf) into H^+^{U{x^),R'^'^). 

We claim that ^2,a = Sjli ^2,j when both sides are viewed as operators from £'{0,, M?'^) to 
P'(C/(x,),M2'^). Indeed, for any / G C^{n,R^'^), since Xj = I on Aq and supp W„ C Aq, 

we have 

w^{x,OfWx,0)ii - XeM\\)) = wUx,^)f{Hx,0){^ - XeMm) ■ {Y.^j{x,0) 

for all x G ?7(x*). Thus 0^2,0/ = Xljli ^2,jf and the claim follows. This completes the proof of 
the lemma. 

Finally, note that 9Tq = OTi^a + ^2,a- Theorem 13.21 follows from Lemma 16.61 and 16.71 

7 Sensitivity analysis of recovering the velocity field from the 
DDtN map 

In this section, we prove Theorem 13.31 on the sensitivity of the inverse problem of recovering the 
velocity field from the DDtN map. We first present a lemma which is a direct consequence of 
Theorem 13.11 and Lemma |3. 31 
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Lemma 7.1. Let c and c be two velocity field in 2l(eo, il, Mq, T), and let f be as in |i^[ ). Then 
there exists 5 > such that if \\Ac — Ac||j:^i[o^3ej/4]xr-5>L2([o,T+€i]xr) — ^> then 

ll^/llL°°(sir,R2'') < C'||/||^i(f^^jg2d) (53) 
for constant C > depending Mq. 

Proof of Theorem 13.31 The proof is divided into the following six steps. 

Step 1. Since c is fold-regular, for each x G il^^, by Theorem 13.21 there exist a neighborhood 
U{x) of X, a smooth cut-off function a £ Cq^{S'LT) and a smoothing operator d\ such that for 
any Uo{x) d U{x) the following estimate holds 

ll/llL2(C/o(a:),R2<*) ^ \\^af\\m(U{x),R2d-) + || ^2,a/ 1| J^i + \\^f\\m{Q,»2d) (54) 

for ah / G L'^{n^^,R^'^). Moreover, both ai2,a and d\ are compact from L2(fi,^,]R2d) 
i/i(C/(x),M2rf). 

We now fix a neighborhood Uo{x) d [/(x) of x for each x. Note that is compact, there 
exists a finite number of points, say xi, X2, ... xm such that C Ujli Uo{xj). Let be the 
operator associated with each point xj. 

Step 2. Denote by H the Hilbert space Hjii ^Ht^(a^j); I^^"*)- We consider the following 
three operators 

Tf = (Otaj^/, DTqj/j •••) ^aj\f /)) 
Tlf = (at2,ai/,ai2,«J,...,%a„/), 
^2/ = (9^ai/, 5^02/' ^^cim/)- 

It is clear that all three operators are bounded from L'^{U^^^) to H. Moreover, Ti and T2 are 
also compact and the following estimate holds 

ll/llL2(f^,R-) < \\Tf\\H + \\Tif\\H + mfWH. (55) 

Step 3. Let £0 C L'^{n,^,R'^'^) be the kernel of T. We claim that £0 C L'^{n,^,R'^'^) is of 
finite dimension. We prove by contradiction. Assume the contrary, then there exists an infinity 
number of orthogonal vectors in £0 C L'^{U^^^,R?'^), say, ei, 62, such that ||ej ||j;^2(q^^ ig2d) = 1 

and Tcj = for all j G N. Since the sequence 

{ej}°Li is bounded in L^{n,o,R'^'^) and the 
operators Ti and T2 are compact, we can find a subsequence, still denoted by {cj}^^, such that 
both the sequences {Tiej}^^^ and {T2ej}'^i are Cauchy in H. By applying Inequality (f55]l to 
the vectors — ej and recall that T{ei — ej) = 0, we conclude that the sequence {ej}^^ is also 
Cauchy in L'^{Q,^^j,R'^'^). This contradicts to the fact that \\ei — ||2,2(q^^ ^2^) > 1 for all i 7^ j. 
This contradiction proves the claim. 

Step 4. Denote by £g" the orthogonal space to £0 in -^^^(^eo' ^^'^)- We claim that 

ll/llL2(n.„M2.) < IIT/IIh for an (56) 

Indeed, assume the contrary, there exists a sequence {fn}'^=i C such that ||/n||L2(nEQ,R2d) = 1 
and ||r/„||j^ < ^ for all n. By the same argument as in Step 3, we can find a subsequence, 
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still denoted by {ej}^-^^, such that both the sequences {Tiej}°2,^ and {T2ej}°2,^ are Cauchy 
in H. By Inequality (i55]l and the fact that ||T/„||/i- < ^, we can conclude that {/n}^i is 
also Cauchy in L'^{n^g,R'^'^). Let /o = lim„^oo/n, then ||T/o||h = lim^^oo \\Tfn\\H = 0. This 
implies that /q G £o- However, note that £,q is closed, as the limit of a sequence of functions in 
SIq, fo must belong to £q . Therefore, we see that /o = 0. But this contradicts to the fact that 
\\fo\\L2{n,^,R2d) = lim„^oo ||/n||L2(!^^^_iR2d) = 1. The claim is proved. 

Step 5. From now on, let / be as in (jl4p . We claim that 



\\Tf\\H 



15 , 



Indeed, for each J^^. , by Lemma [7T1 we have \\'3aj f\\L^(s^r) 

< Apply Ji^ to both 

sides and use the fact that 3a . is bounded from to L^(see [22]), we obtain 



l^aj/llL2(Q^M2d) 



2 



(57) 



Then, 

ll^aj/llHi({/(xj),M2d) 



< 



< 



< 



< 



< 



1 



l^a/ll 



1 
1 

5 
2 



4 
3 

Ci(t7,R2d) 

4 
3 



L2(n,IR2d). 



(by interpolation inequality) 

(by m) 

(by Lemma 16. 6|, I6.7P 

(by interpolation inequality) 

(by interpolation inequality) 



It follows that 



M 



r/ii// = 5^m«,/iiH= 



(C/(a:j),M2d) 



< 



5 



L2(n,IR2d). 



(58) 



to the space L^($7,M'^) by taking 



This finishes the proof of our claim. 

Step 6. Denote by £ the projection of 2q from L'^{Q,R^'^) 
the last three components. Note that the first three components of / are zero, see ()14p . Thus 
the condition V(lnc^ — Inc^) ± £ implies that / G £q. Consequently, Inequality ([56]) holds. 
Combining this with (j58p . we see that 



L2(n,R2d) < 



2 

5 

15 , 



L2(n,i82d). 



Therefore, we must have / = for 



sufficiently small. Finally, note that 



15 , 

■ ell 17 ^ and that both c and c vanishes near the boundary, we conclude that f 



15 , 



implies 



< 
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c = c. This completes the proof of the theorem. 

8 Appendix 

8.1 Linearization of ODE system 

Given the following ODE system: 

y = f{y), y = yo, 

where / G C"'^(M'^, M'^). We consider the pertm'bed the system 

ye = fe{ye), J^e = yo, 

where f^ = f + eg with g G C^{W^). We formally write = y{t) + r{t) = y{t) + e(^(i) + ri(t), 
where i?!)(0) = r(0) = ri(0) = 0. By substituting ye{t) = y{t) + e(t){t) + r\{t) into the perturbed 
system, we can derive that (j) satisfies the following equation: 

m = y^{y{t))-m + 9{y{t)), </)(o) = o. (59) 

By Grownwall's inequality, we can show that \r{t)\ < Ce and < Ce^, where C is a constant 

depending on ||/||c2 + \\g\\c^. 

We solve equation ([59]) as follows. Let A{t) = ^j{y{t)). Let and ^'(t) be the solution 
to the following ODE system: 

m = -mm, m) = id-, 

i!{t) = ^(0) = Id. 

A straightforward calculation shows that ^{t)'^{t) = $(0)^'(0) = Id. Moreover, 

<p{t) = m~^ f Hs)9{y{s))ds. 
Jo 

8.2 Proof of Lemma 14.21 and Lemma 14.41 

Proof of 14. 2i We only show (j27p . since (j28p follows in a similar way. For simplicity, denote 
D = (3ei/4,ti + ei/2) x U{xi). We first show that 

g{t,y,X)=g,it,y,X) + 0{^) in L\D). (60) 

Indeed, by direct calculation, 

git,y,X)-g*it,y,X) = (a(t) - a(ti))e^^^* + a(t)(e*^^ - e^^^*). (61) 



28 



It suffices to show that 

Ri := ||(a(t)-a(ti))e*^^*||i2(B)<^, 

We first estimate Ri. By Lemma 3.1 in [5], we have \a{t)\ « A*. By equation (123]) . we further 
derive that \d{t)\ A*, thus 

a{t)-a{ti)= / a{ti + s{t-ti))ds{t-ti) = 0{Xi)\t-ti\. 



Therefore, 



This proves -^i ^ 

We next estimate i?2- Write f = f* + Sf, then = 0{\{t — ti,y)\^) and hence |1 — e^''^^'^\ < 
X- 0{\{t -ti,y)f). It follows that 

R2< [ \a{t)e'^^*\^-\l-e'^'^\dtdy< [ A^ • A • |(t - ti, y)| V2^(*-*^'?^)^*^(*i)(*-*i'^)' dtdy < 
Jd Jd X 

This completes the proof of (f60]) . 

We now proceed to show (j27p . By direct calculation, 

dg dg^ -^^t ^ . . dn . 

« = ^'^^ ■ 5 - ^A— • 

ay ay ay ay 

df dn dn . ^ 

One can check that ^ — = 0\{t — ti, y)p, then a similar argument as used in the estimate 
of Ri above shows that 



Besides, ([UU|) implies that 



• (9-5*)lli2(m ^ ^• 



Qy • vy - y*7iiL2(D) 
Combining these two estimates together, we conclude that 

ndg dg^ 2 < > 

Similarly, we can show that 
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This completes the proof of ([27]) and hence the lemma. 

Proof of Lemma 14.41 Denote D = (3ei/4, ti + ei/2) x U{xi) again. We first show ([30 
Since x is restricted to V{xi) C T, it suffices to show that 

r(i,y,A) + g(t,y,0,A) = 0(Va) in H\D). 

But this is a direct consequence of Lemma 14.21 and the fact that g = —g*. 
We now prove ([3T]) . By direct calculate 



^{t,x) = ^ia{t)e^'^^''^^)=^Xg■^ 
ou ou ov 

= iXg ■ {^{t) + M{t){x - x{t),u{x)) 

= iXg • {^{h),u{xi)) + iXg • ((^(i), i^{x)) - {C{ti),u{xi))) + iXg ■ {M{t){x - x{t),u{x)). 

Note that in the coordinate x = F{y), 

\{M{t){x-x{t)),u{x))\ = 0{\{t-h,y)\), 
\m,Hx))-{ati),Hxi))\ = 0{\{t-h,y)\). 



It follows that 



\\g-{m,Hx))-m),Hxi)))\\hiD) < - 



\g-{M{t){x-x{t),i^{x))\\l2(^o) < -- 



A' 
1 

A' 



Thus 



Similarly, 



^{t, x) = iXg ■ iah), y{xi)) + 0{y/X). 
-it,x) = iXg~ ■ {r{ti),i^{xi)) + 0{^). 



Finally, using (pOj) and the fact that i^(xi)) = —{(,{ti),h'{xi)), we conclude that ([3T 

holds. This completes the proof of the lemma. 



8.3 An estimate on Gaussian beam interactions 

Lemma 8.1. Assume that Mi and M2 are two symmetric positive definite matrices such that 
< Co < Mi,M2 < ci, and Ni and N2 are two symmetric matrices such that ||A^i||, II-/V2II < C2. 
Let Sx,6$, be two vectors in and A » 1. Then exists C3 > depending only on co,ci and C2 
such that 

,iX-{S^,x)-Xx^ {Mi+i-Ni)x-X{x-Sx)'^ {M2+i-N2){x-5x) i < ^ -C3A(|fap + |3gp) 

A2 

Proof. See Lemma 3.7 in [5]. 
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